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Abstract 

We present a theory of resonant processes in a frozen gas of atoms interacting 
via dipole-dipole potentials that vary as r -3 , where r is the interatomic sepa- 
ration. We supply an exact result for a single atom in a given state interacting 
resonantly with a random gas of atoms in a different state. The time develop- 
ment of the transition process is calculated both on- and off-resonance, and 
the linewidth with respect to detuning is obtained as a function of time t. We 
introduce a random spin Hamiltonian to model a dense system of resonators 
and show how it reduces to the previous model in the limit of a sparse system. 
We derive approximate equations for the average effective spin, and we use 
them to model the behavior seen in the experiments of Anderson et al. and 
Lowell et al. ||||]. The approach to equilibrium is found to be proportional 
to exp(— y/^egt ), where the constant ^y eq is explicitly related to the system's 
parameters. 

PACS numbers: 34.10. +x, 34.90.+q 
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I. INTRODUCTION 



Frozen gases are a new and in some ways ideal laboratory to test our understanding of 
quantum theory in a complex system. With present technology one can manipulate and 
detect electronic processes with an extraordinary selectivity and precision; furthermore, the 
translational temperature of the gas can be lowered to the point where it can be ignored 
when discussing electronic processes. Frozen Rydberg gases have the added advantages 
that their states are well understood and many processes occur on a microsecond time 
scale, easily allowing for time-resolved spectroscopy. Pioneering experiments on resonant 
processes in these gases have been carried out by Anderson et al. |TJ and Mourachko et al. 
@. The present work was motivated by the desire to understand those experiments and the 
subsequent work of Lowell et al. |3|,|3J highlighting the dynamic aspects of these resonant, 
many particle systems. The measurements are of mixtures of S5 Rb atoms initially prepared 
in the 23s and 33s states (henceforth to be called the s and s' states, respectively). There 
are initially N atoms in the s state and N' atoms in the s' state. The transition ss' — > pp' 
is monitored, where p refers to the 24p state and p' refers to the 34p state, at and near 
resonance, where e' s + e s = e' p + e p . Some of the experimental features we wish to understand 
are the rapid rise followed by slow approach to saturation of the signal and the width of the 
lineshape (in the detuning A = e' p + e p — e' s — e s ) as a function of time. 

Each individual interaction leads to a (coherent) oscillatory behavior, but as we will see 
below, averaging the ss' — > pp' interaction over the random positions of the atoms greatly 
smooths out the signal. The effective "incoherence" brought about by the sp — > ps process 
(the "walking away" discussed by Mourachko et al. @) completes the smoothing out of the 
on-resonance signal but has less effect on the off-resonance signal. (The s'p' — > p's' process 
is also present, but is smaller than sp — > ps by a factor of 16 in this case and can therefore 
be neglected to the first approximation.) 

It was surmised by Anderson et al. [|J that the linewidth must be of the order of the 
average interaction energy. This interaction splits the ss'-pp' degeneracy and the subsequent 
migration of the p state to other atoms broadens the energy of this "elementary excitation" 
into a band, as appropriate to an amorphous solid. Cluster calculations were performed 
to illustrate this band formation ||, but averaging over atom positions was done by order- 
of- magnitude arguments only. Mourachko et al. |2j pointed out, in connection with their 
experiments on a different system, that the "walking away" of the p excitation from its 
original location should be regarded as a diffusion process, and dynamical equations for the 
resonance transition in the presence of this diffusion were written down ||. 

The present work builds on these earlier insights, but goes considerably further in the 
process of accurately averaging over atomic positions. It turns out that this averaging 
(effectively, a phase averaging) itself produces an exp(— y/jegt) dependence at large times, 
characteristic of a diffusive process. A report of the newer experimental data by Lowell 
et al. in comparison with the results of the present theory is being submitted separately 
In this paper we give a full account of the general theory, and we discuss in detail in 
Section |I| a particular case that cleanly shows the effects of randomness and phase-averaging. 
In this simple case, the averaging process can be carried out exactly. For the practically 
important r~ 3 interaction potential, one can proceed in a mathematically elegant way and 
obtain closed formulas that can be evaluated by a computer package, or in some cases by 
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hand. It appears to be a fortunate coincidence that the r~ 3 potential accurately describes 
the dipole-dipole interactions that are so prevalent in nature. Some aspects of this approach 
are easily extended to more general interactions; however, the r~ 3 dependence has special 
properties which, when coupled with the assumption of a random distribution of atoms, lead 
to relatively simple results. 

Although this work was stimulated by experiments on Rydberg atoms, it is generally 
applicable to resonances induced by dipole-dipole interactions. These may be common in 
highly excited gases and in molecular systems. Frozen, resonant gases may even exist in 
interstellar clouds. 

We develop the theory in three stages. First, in Section [Tl], we discuss a single s' atom 
interacting resonantly with a surrounding gas of s atoms, without any possibility for the s 
excitation to "walk away" through an sp — > ps process. When the result is averaged over the 
positions of the s atoms, it corresponds to a "sparse" system of s' atoms, i.e. a system where 
N' <C N. We do not know of an existing experiment to which this treatment applies, so this 
section can be viewed as a theoretical prediction of the outcome of a possible experiment. 
Our purpose here is also one of exposition: this example allows us to introduce in the simplest 
context some of the mathematical techniques used throughout the paper. In Section |T|, 
we complicate things somewhat by allowing the s state to decay away exponentially (or 
more generally, but explicit formulas are given for the exponential decay). If the s state 
decays diffusively, our treatment here is related to that of Mourachko et al. @, except for 
the introduction of mathematically exact phase averaging. Here, too, the treatment applies 
only to the case where N' <C N. Finally, in Section [TV] we introduce an effective spin 
Hamiltonian which fully models the atomic system for all values of N and N', as well as for 
any strength of the resonant process (ss' — > pp') and of the mixing processes (sp — > ps and 
s'p' — > p's'). This system is relevant to the experiments of Anderson et al. |[J and Lowell et 
al. 0,3- We show its equivalence to the previous results in the sparse limit (N' <C N) and 
discuss briefly the approximate solution in the non-sparse case. The results in this section 
are summarized elsewhere ||, but the derivation is presented here for the first time, as are 
additional formulas. The limitations of the present theory and some of the many possible 
extensions are discussed briefly in Section |V|. 



II. SPARSE SS' -» PP' PROCESSES 

We consider one atom, initially in the state s', in interaction with a gas of s atoms 
through an ss' — > pp' process. Like Mourachko et al. ||, we describe the system by the 
equations: 

i d = A a + Vk c k , (la) 

k 

i c k = V k a . (lb) 

Here ao(t) is the amplitude of the state in which the atom at the origin is in state s' and all 
other atoms are in state s, while c k (t) (with k running from 1 to N) is the amplitude of the 
state in which the atom at the origin is in state p' and the atom at is in state p, while all 
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the others remain in state s. Vk is the interaction potential and A = e' p + e p — e' s — e s is the 
detuning from resonance. We will assume that is of the dipole-dipole form 

V t = 4 (2) 



r 



The atoms are sufficiently cold that during the time scale of interest they move only a very 
small fraction of their separation, and therefore Vk can be taken to be independent of time. 
As the temperature increases one approaches the opposite limit, where binary collisions 
control the resonant process |§. There must of course be a gradual transition, which is 
beyond the scope of this paper. 

The set of differential equations (1) can easily be solved with the initial condition a^it = 
0) = 1 and Ck(t = 0) = for all k, yielding 



2i V k e~ iAt / 2 sin(v / A 2 + 4V 2 t/2) 
^ (*) = 7==^ — , (3) 



VA 2 + 4V 2 

with 



N 

v 2 = E^4 2 - (4) 



fc=l 



Eq. (|Tb|) then gives a (t). The experiment monitors an optical transition from the p' state 
(or, equivalently, from one of the p states), which is proportional to 



S(t) = l-\a Q (t)\ 2 = J2\ck(t)f 

k 

4V 2 sin 2 ( v / A 2 + 4V 2 t/2 
= A 2 + 4V 2 



(5) 



We see that the signal S (t) exhibits Rabi oscillations. 

We apply this result to a system of N' atoms, initially in state s', randomly dispersed 
among a much larger number, N, of s atoms. In this sparse limit, the signal is proportional 
to N' times the sample average of S(t), and the sample average is equivalent to an ensemble 
average over the atomic positions that are hidden in V 2 (see Eq. (0) below). In general, 
averaging is more easily done on the Laplace transform of the signal, which in this case is 



_ roc 

S(a) = / e- at S(t) dt 
Jo 



1 2V 



2 



a a 2 + A 2 + 4V 2 
1 / a 2 + A 2 



2a y 1 a 2 + A 2 + 4V 2 J ^ 

Later in this paper, more complex models are solved directly by Laplace (or Fourier) trans- 
forms; thus Eq. (|J) is useful for comparison with these more general results and approxima- 
tions. 



4 



A. Averaging over atom positions 



To compute ensemble averages, we use the following result, valid in the limit N ^> 1, for 
a random distribution of the variables r^: 



N 



k=l 



^ / d 3 rexp -/3V 2 (r) 



N 



1 



exp 



1 

ft 



d 3 r 



N 



d 3 r 



-f3V 2 (r) 
_ e -PV 2 (r) 



N 



(7) 



where Q is the volume of the gas in the trap. In particular, for the dipolar interaction 
V(r) = fifi'/r 3 , we have 



d 3 r 



I _ e -/J(wO ! V 



leading to 



where 



\ I ( 



4^3/2 



-vWf3 



4tt 3 / 2 N , 

— n ^ 



(8) 



(9) 



(10) 



For the typical densities in the experiments of Anderson et al. 0] and Lowell et al. [|3],U , v is 
on the order of MHz. (In Appendix |A], we repeat the arguments above with an interaction 
that includes the angular dependence due to the relative positions of the dipole moments and 
assumes that all of the dipole moments point in the same direction.) Since v has the units 
of an energy, it must be of order fifi'N/Q on dimensional grounds. It is still remarkable, 
however, that Eq. ( |10"D gives simply and exactly the quantity v that will enter in all the 
averaged quantities in this section. 

Using Eq. (|9]), we can evaluate the average of any function F(V 2 ) provided that it can 
be represented as 



F(V 2 



d(3 F{(3) e-^ 2 



Further, we can use the following identity 



7T JO 



dy exp 



v_ 

' 4 



v 2 (3' 



(11) 



(12) 



so that 
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F(V 2 ) av = / d/3F(/3)e-V/3 
Jo 

= -4= / dye"" 2 / 4 / df3F(f3)e- v2/3/y2 

\/7T Jo Jo 



■^r*^(5)- (i3) 

This last set of relations implies that we can average a function F(V 2 ) over the positions 
of the interacting dipoles (atoms in our case) by replacing V 2 with v 2 /y 2 and integrating 
over the kernel exp (— y 2 /A) /y/n. This trick is not always useful, as it can lead to highly 
oscillatory integrals. It is, however, convenient in determining the lineshape in the saturation 
(t — ► oo) limit (as shown below). 



B. The averaged signal and lineshape 

Starting from Eq. (H) in the form 

D -l3(a 2 +A' 2 +4V 2 ) 



S(a) 



a 2 + A 2 



2a 



2a Jo 



dp e 



(14) 



and using Eq. (0), we obtain 



2a 2a Jo 



7T v 



2 aA 



exp 



v 

A* 



erfc 



(15) 



where A 2 = a 2 + A 2 . Expanding in v/A and performing the inverse Laplace transform leads 
to 



S(t) c 



-vt) n 



2 rf Sr(n + 2)r(f) 



1-^2 



n + 1 n + 2 n + 3 A 2 t 2 



2 ' 2 



(16) 



where n F m (ai, . . . , a n ; b%, . . . , b m ; z) is a generalized hypergeometric function. It is known 
that iF 2 is related to the integral of a power multiplied by a Bessel function. Expanding 
iF 2 yields the double series 



oo oo 



s(t) av = ^EE 



n=0 m=0 



4 m n\m\T 



+2m+2 



(n + 2m + 1) 



(17) 



which can also be obtained directly by expanding Eq. (|15D in l/oi and then taking the inverse 
Laplace transform. The double series can be rearranged to obtain the following expansion 
in the detuning: 



2 Jo ^/p ^ ml 



2v 



J2m+1 V 8^pt . 



(18) 
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The first term in this expansion, the on-resonance signal, can be written as 



S(t) c 



A=0 




(19) 



which can also be obtained by setting A equal to zero in Eq. (|T6| ) and summing the resulting 
series. Plots of S(t) av for several values of A are shown in Fig. |l|. Notice that the initial 
slope of the signal is independent of A. This feature translates into resonance widths that 
vary as 1/t for small t - the so-called transform broadening discussed by Thomson et al. 
With v on the order of MHz, this initial rise occurs in a fraction of a [is. Another point of 
interest is that the averaging is much more effective at smoothing out the oscillations for 
A = than it is for A^O. 

We can construct a measure of the width from the first two terms in the A-expansion, 
Eq. ©.If 



and we define 



S(t) av = S (t) + A 2 ^ (t) + O ( A 



(20) 



w 



S (t) 



(21) 



then 2w would be the FWHM if the lineshape were Lorentzian. The quantity w/v is plotted 
in Fig. [2]. For small vt, we have explicitly w = y/l2/t. While the 1/t behavior follows 
from the transform broadening argument, or even simply from dimensional analysis, the 
coefficient VT2 is a prediction of the detailed theory. 

The lineshape at saturation (t — > oo) can be obtained using the relation in Eq. (|13|). 
which yields 



1 f°° 2 

S(t) av = —= / dy e y 

\/7T JO 



l-cos U/A 2 + ^ t 



1 I A 2 9 



For large times the cosine term above averages to zero, giving 

,2 \ 



(22) 



S(t 



oo 



2 A 6XP U erfC U 



(23) 



which can also be extracted from the small a behavior of Eq. flToQ. This lineshape is plotted 
in Fig. ^. The FWHM is approximately 4.6?;. Also plotted is a Lorentzian lineshape with 
the same height and FWHM. Notice that this lineshape is sharper than the Lorentzian for 
small A and falls off more slowly for large A. 
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III. SPARSE SS' -> PP' IN A BATH OF SP -> PS 1 

We again consider one atom, initially in state s', in interaction with a gas of s atoms via 
the ss' — > pp' process, but we now allow for sp — > ps processes by modifying Eqs. (1) to 

i d = A a + J2 V k c fc> ( 24a ) 

k 

i c k = V k a + J2 U k ici . (24b) 
i 

where U^i = jU 2 /rjy and the other symbols have the same meaning as in Eqs. (1). Although 
Eqs. (24) are exact, they are difficult to solve exactly. We therefore introduce the equation 

% c k = V k a - i'jCk , (25) 



in place of Eq. ([245), where 7 is an effective inverse lifetime. In a crude way, the — ryc/c 



term in Eq. (|25|) describes the "walking away" of the p states from the neighborhood of the 
resonating atom at the origin. One can see from Eq. ( [24 b| ) that the model of Section |I| is 
recovered in the limit \x <C //. 

When Eq. flTED of the last section is replaced with (fffi), the solution of the new system 
can be obtained by replacing A with A + 27 in Eq. (^) and in the corresponding equation 
for do(t). Then 



/V 



£(*) = £M*)f 
fc=i 

2V 2 e" 7 * [cosh (yt) - cos (xt)] 
x z + y 2 

where x and y are real numbers that satisfy 

(x + iyf = (A + i 7 ) 2 + 4V 2 . (27) 

At resonance, we have 

= I (^^^ <2V - (28) 
V ' y; \ (0, V7 2 - 4V 2 ) if 7 > 2V. 

Note that L (t) decays more slowly in the latter case. Hence the large-t limit of L (t) av will 
be dominated by small V. This remains true off-resonance as well. The small-V expansion 
for y is 

Thus, we can conclude that for the leading behavior at large times, 

2V 2 7 



lim L (t) oc exp 

t— >oo 



A 2 + 7 : 



-t 



(30) 
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Using Eq. (Q), we obtain an asymptotic behavior of the average 



Km L {t) av oc exp (-^7^) , (31) 

where 

7e 9 = ^p^- (32) 

Let us now examine the full time dependence. We begin with the Laplace transform of 
L{t) which is 

f, x_ 2sV 2 
W s 4 + (x 2 - y 2 ) s 2 - x V 

where s = a + 7. Substituting x and y from Eq. fl27|) gives 

2sV 2 

= S 4 + ( A 2_ 7 2 +4V 2) S 2_ A 2 7 2- ( 34 ) 

The signal is still proportional to N'S(t), with S(t) = 1 — |a (t)| 2 . It is no longer true 
that 1 — I do 1 2 = J2k \ c k\ 2 , because J2k \ c k\ 2 is n °t the total probability of finding ap state, but 
only the probability of finding a p state that has not "walked away" . However, the relation 



|(|ao| 2 + El^lj=-27El^l 



(35) 



with the usual initial conditions, enables us to obtain 

S(t)=L(t) + 2j [ L(t')dt', (36) 
Jo 



and 



S(a) = [ 1 + — ) L{a). (37) 



a 



For 7 > 0, S(t) rises from to 1 (or more accurately to (1+ N'/N)' 1 , as seen in the next 
section), while for 7 = it saturates at 2V 2 / (A 2 + 4V 2 ). What happens for small 7 is a 
quick rise (on the average time scale of v ~ 1 ) to 2V 2 / (A 2 + 4V 2 ), followed by a slow rise to 
1. 

Averaging these equations is done as in Section ||. For instance, averaging L (a) yields 



7T v 1 v 2 \ „ 1 v 



L ( a )«« = 2sB 6XP \lp ) 61(0 \b) ' < 38) 
where 

5 2 = (V + A 2 )(V-7 2 )/ S 2 , (39) 
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which is the analog of Eq. (|15D . The inverse Laplace transform of this expression is rather 
complicated, but it simplifies immensely at resonance, becoming 

* w- - f «f (f)" /2 " 7 'i ,<iTT " /2/ "' 2 M ■ (40) 

This result can also be expressed in terms of generalized hypergeometrics, and in that form 
appears similar to Eq. (|16|). One can use the asymptotic behavior of the Bessel function 
in Eq. (KJ) to confirm the large- 1 dependence argued earlier in Eq. (|3"T|). An expansion 
in inverse powers of a, leading to an expression similar to Eq. (|17D , is easily obtained by 
computer. Fig. § shows the average signal S (i) for 7 = v , and A/v having the values 0, 1, 
and 2. One can see that the approach to one becomes slower as the detuning is increased, as 
predicted by Eqs. (|3lf) and fl32|) . This behavior is the same as that found by Mourachko et 
al. [0 under the assumption that the time development of the is governed by a diffusion 
equation. We see that in the present approach it arises simply from phase averaging. The 
result for the on- resonance signal is shown in Fig. ^|. As one would expect, for small values 
of 7 the small-t features of Fig. |^ resemble those in Fig. [l| of Section U, where the sp — > ps 
process is neglected entirely. 



IV. INTERACTING RESONANT PROCESS 

We now consider N atoms at positions r& which are initially in state s and N' atoms at 
positions ry which are initially in state s' . We represent the s and p states at with the 
down and up states of an effective spin a>, and similarly represent the s' and p' states at 
with a spin ay. The Hamiltonian is then 



J2 [e. + (e p - e s )a+a k ] + [ e - 



{e' p - e' s )a£,a k , 

k k' 



+ J2 V kk> [vtet' + o k °k\ + J2 U ki°t°i > ( 41 ) 

kk' ~ k,l^k 



where 



/ 2 



Vkk' = 1 ^3 and U M = -3. (42) 

\ r k — r k'\ \ r k — r l\ 

For the experiments under consideration \i w 4/i', and the sp-ps coupling /i 2 , which leads to 
"spin diffusion" , is larger than the ss'-pp' coupling ////, which is responsible for the resonant 
energy transfer. 

The strategy is to write down the evolution equations for the spin variables, take the 
expectation values over the initial state \i) (which consists of all spins down), discard the 
expectation values of quantities that fluctuate incoherently, solve the resulting equations, 
and finally average over the atomic coordinates and ry. The experimental signal is 
proportional to the number of p' states, J2y n k'i where 

n W = ~ (1 + K,)) ■ (43) 
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It should be noted that \i) is not the ground state, nor does the system relax to the ground 
state, because it is essentially decoupled from other degrees of freedom for the duration of 
the experiment. 

We start with basic equations of motion such as 

ij°t> = A 4 + °t> E Vk>k<Tk, (44) 

k 

= oi E V kk'Ok> + o'l E Uuvt, (45) 

al k' l+k 

l J t a l> = 2 E V*k °% ~ °k°k) , (46) 

where A is the detuning from resonance |J, and obtain 

1 d 2 

2^2°"^' = ~ ^ V k'kV k i'Cr z k (o k ,aj + a%a\, 

kV 



- E Vk'kVwoi' (°k °i + °t°, 



ki 



- E V k'kUki(? z k [a^af + a k> a i 



-AE^K^ + ^j- (47) 
k 

We reduce the expectation value of this equation as follows. In the first line, we assume that 
°% {®k' a v + °k' v\) * s negligible unless k' = I', in which case it reduces to (a k ) . Similarly, 

a k , (a^af + G k a lj) ^ s replaced by (a k ,) 5 k i. The third line contains only incoherent terms, 
which we treat by assuming that they lead to a Lorentzian broadening by letting 

< E UutT? - ilo+ (48) 

l^k 

when computing these terms. With this crude approximation and Eq. fl4"6|), we obtain 



E Vk'kUu (a* h (aW + w)) = lj t Wl ■ (49) 



ki 

Collecting results, we obtain 

K>) =-E^'^'(K) + K')) 

-lj t {a z k ,)-Ag k ,, (50) 

where 

9k'(t) =J2Vk>k(a£,a£ + a k ,a k ). (51) 
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Using arguments similar to those used to derive Eq. (|50|), we obtain the approximate equation 

$-A£- 7fa ,-», (52) 

where g = — Ylk' 9k' an d the 7$ term has been put in by hand to preserve conservation. 
These equations imply 

d 9 A dn ' (ro^ 

Tt =A lF> (53) 

where 

™' = ^E n ^ ( 54 ) 
iv k' 

with the initial condition gk'{0) = 0. 

At this point, it is more convenient to work with the p' occupation n^, defined in Eq. (fl3|), 
and with the corresponding p occupation, n&. For comparison with experiment we will need 
the sample average, which corresponds to S (t) of Sections [TJ and fTJ. Taking Laplace 
transforms of Eqs. (|50D and (0), we obtain 



2 




„ E " 2 E - A 2 -f-n', (55) 

fe A: a + 7 

where the tildes denote Laplace transforms and n'{a) is the sample average of hy(ot). 

We can obtain a similar equation for n^, and we are then left with coupled equations 
which we are not able to solve. However, if N pa N' we can simplify the problem by assuming 
that the strong sp-ps coupling acts to randomize the spins, while the resonant ss'-pp' 
processes are entirely coherent. Accordingly, in Eq. ( ]5"5] ) we replace with its sample 
average n and obtain an equation for alone by using the exact relation Nn = N'n', 
which follows from the fact that p and p' states are created in pairs by the Vkk' term in 
Eq. PI]). We solve for and average over a random distribution of atoms, replacing the 
sample average (1/N') J2k ,n k' with an ensemble average, and obtain 



n' 



[l/a) (l - C 2 F a . 



1 _L N' _ (N> rt-2 _ jA2_\ fp 
L ~T~ N V a+7/ a 

where 



(56) 



A 2 

C 2 = « 2 + « 7 + ^^, (57) 
a + 7 



and F av is the ensemble average of 
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For N = N', Eq. (M) simplifies to 



n = : .a^/ . ( 59 ) 



2-(C 2 -^F, 



ar 



F av is computed in the same way as the averages of S in Section |T| and L in Section [HI 
The analog of Eq. (|15]) is 



7T f / f 2 



C^U = l-^exp^je rf c^j, (60) 
with v given by Eq. ([TDD and C by Eq. fl57|). One way to compute n'(i) is to expand h'(a), 



Eq. (^), in inverse powers of a. This gives the analog of Eq. (17D, but we have not found a 
simple expression for the coefficients. 

Fig. |] shows the signal as calculated in Eq. (|59|) for 7 = v , and A/v — 0, 1, and 2. We see 
that the on-resonance signal is devoid of any oscillations. Off-resonance, the oscillations are 
not completely washed out by the phase averaging, and some nonmonotonicity is evident. 
Consequently, the dip in the width seen in Fig. ^| is also present in this calculation, provided 
7 is not too large. 



Fig. |6] differs from the corresponding Fig. ^ of Section |T| in one important respect: 



the on-resonance signal saturates to 1/2 rather than to 1. For general values of N and 
N', Eq. (|56|) predicts that as t — > 00 (i.e., as a — > 0) N'n' approaches NN' / (N + N'), as 
expected from simple kinetics for the two-body reaction ss' — > pp'. 

As in the previous sections, the energy scale is set by v, which according to Eq. (|T0|) has 
the value 5.72 fifi'N/Q. As discussed by Anderson et al. this quantity is of the correct 
order of magnitude to account for the data. 



A. Reduction to the sparse limit 

We now show explicitly how the general spin- variable formalism relates to the coefficients 
ao and c& of Sections |TT] and [TTI| when there is only a single atom that can be in the s', p' 
pair of states (i.e., when only a single spin variable of the type ov is present.) The key to 
the correspondence is that in this case the initial state \i) (which consists of all spins down) 
evolves to 

I*>=(oo(*)+E c «(*)4<jl0 (61) 

at time t, where the spin variables now denote time-independent (Schrodinger) operators. 
With the understanding that = 0, we have: 

K,) = (t|^|t) = l-2|a (t)i 2 , (62) 
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(a z k ) = 2\c k (t)\ 2 -l, 



(63) 



(a+a+) = a (t)c{(t). (64) 

It is then clear that the expectation value of Eq. ([5]) is equivalent to the following, which 
is itself a consequence of Eqs. (1): 

■j t KWI 2 = V k ( fl o4 - a l C k) , (65) 

where Vk is short for Vk'k- Similarly, the equations for (cr|) and (o'k' a k} are equivalent to 
the equations for |cfc(t)| 2 and a (£)cfc(£)t that follow from Eqs. (|24b| ) . 

When the Lorentzian approximation is made in Eq. (|47|) or its analog for <j|, it is possible 
to obtain a closed set of equations for (a%) and (of,) that give the results reported in 
Sections [TI| and [TIT] for a sparse system. As indicated by the Laplace transforms in Eqs. ( |34]) 
and (|37|) , the equation for L{t) = \ (1 + J2k ( a k)) mus t be of the fourth order, and that for 
S(t) = \ (1 - Efc (<**)) of fifth order. 

If Eq. (|56|) were exact, it would correctly give the signal in the limit N' <C N, and n'(£) 
would be equal to S(t) av of Section [TTJ if the values for 7 used in the two models were assumed 
to be equal. In addition, n'(t) would reduce to the exact S(t) av of Section [Ij if 7 were taken 
to be zero. Since the result for n'(t) is based on approximations that are certainly not valid 
for N' <C N, it is already comforting that, in that limit, it shows a general resemblance to 
the correct behavior. We are currently working on better approximations to apply to the 
spin-variable formalism. 



V. DISCUSSION 



The system we are studying can exhibit a wide range of behaviors, depending on the 
values of the parameters fi/fi' and N/N'. At the same time, there are features that persist 
(at least qualitatively) throughout the range of parameters. 

In Section 0, we determine the exact behavior of a frozen system consisting of a single 
s' state interacting resonantly with a sea of s states. That is to say, we treat exactly the 
limit /i fi' and iV' C iV. A single ss' pair produces, of course, a signal with a sinusoidal 
behavior, and the same is true, less obviously, for a single s' interacting with a sea of s 
that are not mutually interacting. The phase interference brought in by averaging over 
configurations causes the signal from a sparse system of s' to take on the form shown in 
Fig. [TJ. Because the averaging does not completely wash out the oscillations at small times, 
especially for nonzero values of the detuning, the resonance linewidth has the nonmonotonic 
behavior seen in Fig. ||]. The initial rise of the averaged signal is linear in time, even though 
the unaveraged signal is quadratic, and is given by y/nvt/2, with v defined in Eq. (|T0|). The 
linewidth at large times is approximately 4.6 v ||10|| . The fact that v is sizably larger than 
the simple estimate (N/fl) reflects quantitatively the fact that "close pairs" of atoms 
are more heavily weighted in the average. 
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In Section [HI| we try to allow for any ratio fi/fi', but we still keep N' <C N: physically, 
we consider a sparse distribution of (initially) s' states in a gas where sp-ps flips take 
place at a non-negligible rate. The first question we should answer is: after the ss' — pp' 
transition, does the p state "walk away" to infinity or remain localized? In this paper, we 
have assumed that localization does not take place, or, if it does, is important only in the 
very sparse limit where the localization length is smaller than (il/N) 1 ^ 3 . Leaving much 
room for future work, we have simply introduced a lifetime 7 _1 for the p state to "stick 
around". The spectral function for the propagating sp-ps flip is then simply a Lorentzian. 
We have obtained tentative estimates of 7 when it is entirely due to the sp-ps flips, and 
we plan to report these results in the near future. There can also be other contributions to 
7, arising for instance from translational motion of the atoms. (The gas can be treated as 
frozen to a good approximation when the typical collision time r exceeds v~ l , and certainly 
t -1 contributes to 7.) In this paper, we in effect add a damping term to the equations of 
Section ||, making the problem soluble again. The effect of 7 is to further smooth out and 
wash away the oscillatory behavior of the unaveraged signal. For large 7 all traces of the 
oscillations are wiped out and the signal saturates according to exp (—y/jegt), where ^ eq is 
given by Eq. fl32|). The saturated linewidth is infinite, because the s' state eventually decays 
for all values of A. However, the signal at finite t varies significantly with A, indicating 
that the linewidth has a minimum at intermediate times. For small t the signal is always 
independent of A, consistent with a linewidth proportional to t™ 1 , as seen already in Fig. |2|. 

In Section [TV] we introduce an effective spin Hamiltonian to describe the resonating frozen 
gas for all values of the parameters and N/N', and we discuss in particular the case iV « 
N', which is relevant to the experiments carried out so far. The only simplification we make 
in this spin Hamiltonian is to neglect the potential responsible for the s'p'-p's' interaction. 
This coupling can easily be included, although the resulting equations of motion contain 
more terms and would be more difficult to handle; it is in fact small for the experiments of 
Anderson et al. [|T]|| . Other resonant systems can be described by similar spin Hamiltonians. 
The experiment of Mourachko et al. , for example, consists of a pp — ss' resonance in the 
presence of sp-ps and s'p — ps' flips. A spin-one Hamiltonian therefore describes this system, 
with the —1, 0, and 1 spin states corresponding to the s, p, and s' atomic states, respectively. 
In this case, there is only one set of spins and all spins are initially in the state. The spin 
model we have discussed in this paper consists of two sets of spins, which complicates matters 
but at the same time allows one to consider the limit in which one set is sparse. 

There is practically no end to the variety of phenomena that this type of spin Hamiltonian 
can describe, combining features from many branches of physics. Atomic physics, random 
systems, spin systems and magnetic resonance, and many-body theory are all represented. 
The equilibrium behavior poses interesting problems of statistical mechanics, but we are 
interested in the quantum dynamics, starting from a prepared state. We have obtained 
and solved approximate dynamical equations for the case A" ~ N', using the Lorentzian 
approximation as in Section |T| to model the sp-ps interaction. Much more work can be 
done and is now in progress, and some of the results may change when better approximations 
are introduced. In particular, in a better theory the linewidth is likely to depend on the 
sp-ps interaction as well as on the ss'-pp'. Nevertheless, the qualitative behavior of the 
signal is common to all approximations we have tested so far and matches what is seen in 
the experiments that initially motivated this work. The important features that we already 
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discussed are all there: the sharp initial rise, the slow saturation, the dependence on A that 
is less pronounced at small times and corresponds to a time-dependent linewidth. 
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APPENDIX A: ANGULAR AVERAGING 



In Section |TT| we averaged the quantity e _/3y2 , where V 2 = J2k=i an d Vk — A i / i '/ r i- m 
this appendix we include the angular dependence of the dipole-dipole interaction assuming 
that all of the dipoles point in the same direction, i.e. we replace Vj. with 



The analog of Eq. 



V{r,6) 



(3 cos 2 e- 1) . 



now has a nontrivial angular part; it becomes 

i-tt/2 



2tt 



r dr 



sin 9d9 



n/2 



1 



-/3V 2 (r,e) 



16tt 3 / 2 
9^ 



HH 1 V p. 



(1) 



(2) 



The additional angular dependence changes v by a factor 4/3V3. 
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FIG. 1. The averaged signal, S(t) av , for A/v = 0, 1, and 2 (the upper, middle, and lower 
curves, respectively) calculated from Eq. ( |T7| ) for the sparse limit. The ss'- pp' process is averaged; 
the sp-ps process neglected. The Rabi oscillations are not completely washed out by the averaging 
process and are more pronounced off-resonance. 
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FIG. 2. The "width" versus time, as given by Eq. (^H), shows a dip which is associated with 
the stronger and higher frequency oscillations seen at larger values of A as shown in Fig. |l|. 
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FIG. 3. The saturation (t — > oo) lineshape, Eq. (|23|). It is not a Lorentzian but a Lorentzian 
averaged as in Eq. (|2^). The dashed line is a Lorentzian with the same height and FWHM, drawn 
for comparison. 
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FIG. 4. The averaged signal, S(t) av , for 7 = v, and A/v = 0, 1, and 2 (the upper, middle, and 
lower curves, respectively) calculated from S (a) av = (l + {&) a v an d Eq. (|38[). The approach 

to saturation goes like exp (—y/jeqi), with j eq given by Eq. (|32|). 
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FIG. 5. S(t) av , as in Fig. ||, calculated at resonance (A = 0), for inverse lifetimes j/v = 1, 
0.1, and 0.02 (the upper, middle, and lower curves, respectively). For the largest value of 7 the 
oscillation is completely gone. 
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FIG. 6. The averaged signal, n'(t), for N = N' , ^ = v, and A/v = 0, 1, and 2 (the upper, 
middle, and lower curves, respectively), as calculated from Eq. d59|). While the oscillations are 
"damped out" for A = 0, they remain for A/v = 1 and 2. 
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